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Abstract. Bragg scattering in sinusoidal PT-symmetric complex crystals of finite 
thickness is theoretically investigated by the derivation of exact analytical expressions 
for reflection and transmission coefficients in terms of modified Bessel functions of first 
kind. The analytical results indicate that unidirectional invisibility, recently predicted 
for such crystals by coupled- mode theory [Z. Lin et a/., Phys. Rev. Lett. 106, 
213901 (2011)], breaks down for crystals containing a large number of unit cells. In 
particular, for a given modulation depth in a shallow sinusoidal potential, three regimes 
are encountered as the crystal thickness is increased. At short lengths the crystal 
is refiectionless and invisible when probed from one side (unidirectional invisibility), 
whereas at intermediate lengths the crystal remains refiectionless but not invisible; for 
longer crystals both unidirectional refiectionless and invisibility properties are broken. 
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1. Introduction 

Physical phenomena described by reduced or effective non-Hermitian Hamiltonians are 
often encountered in a wide class of quantum or classical systems, for example in nuclear 
or condensed-matter physics of open systems ^ E] or in optical systems in presence of 
optical gain or losses [4J. Among non-Hermitian Hamiltonians, great interest has been 
devoted in the past two decades to study the properties of parity-time (VT) invariant 
Hamiltonians, which possess a real- valued energy spectrum below a symmetry-breaking 
point in spite of non-Hermiticity. Such a class of non-Hermitian Hamiltonians has been 
originally introduced by Carl Bender in the framework of non-Hermitian extensions 
of quantum mechanics and quantum field theories [51 EJ [7j, and found recently an 
increasing interest since the proposal of physical systems described by T^T-symmetric 
Hamiltonians, including optical P[9l[ini[IIl[l2l[l3l[H[l5l[ini[I3 electronic 
[20] systems. 

Complex periodic potentials [ai[22l[23l[2l[25l[26l[23 reahze a kind 

of synthetic complex crystals, which show rather unusual scattering and transport 
properties as compared to ordinary crystals. Complex crystals have been investigated 
in diflFerent areas of physics, ranging from matter waves [25l [26l [271 EHl ES] to optics 
[To] . In optics, a complex crystal with VT invariance is realized by introduction 
of index and balanced gain/loss modulations in a dielectric medium [LQ\. Complex 
crystals can be also realized in atom optics experiments exploiting the interaction of near 
resonant light with an open two-level system. Scattering of matter waves from purely 
absorbing optical lattices was reported in a few earlier experiments Refs. [251 [261129 ]. As 
noticed by Berry [31j, although in such experiments on matter waves VT symmetry is 
not strictly realized and the complex potentials are purely absorptive, the analysis is 
essentially the same, since the mean loss simply represents an overall exponential decay 
of the wave. From the theoretical side, Bragg scattering, diffraction and transport 
properties have been extensively investigated for sinusoidal PT-symmetric complex 
crystals [TOl [211 [22l [32l [33l [34l [35l [36], revealing some interesting properties such as 
violation of the Friedel's law of Bragg scattering [26l [271 132], double refraction and 
nonreciprocal diffraction [lOj, and unidirectional Bloch oscillations In particular, 
in a recent work [36] it was predicted that a sinusoidal T^T-symmetric sinusoidal 
crystal of finite length near the spontaneous T^T-symmetry breaking point can act as a 
unidirectional invisible medium, i.e. the crystal is almost refiectionless when probed from 
one side, and transmission occurs as if the crystal were absent. Such an unidirectional 
invisibility of T^T-symmetric Bragg scatters near the symmetry-breaking point was 
previously predicted to occur in Ref.[37j for waveguide Bragg gratings which combine 
matched periodic modulations of refractive index and loss/gain yielding asymmetrical 
mode coupling. In these previous studies [36l [37], invisibility was explained on the 
basis of a coupled-mode theory describing Bragg scattering and coupling of counter- 
propagating waves in the crystal, which is rather common in the optical context 
[381 139] . Such an analysis predicts that, for a shallow grating near the VT symmetry 
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breaking point, the sinusoidal crystal appears to be invisible when probed from one side 
independently of the crystal length. 

In this work we re-consider the scattering properties of the sinusoidal T^T-symmetric 
potential and derive exact analytical expressions for reflection and transmission 
coeflicients. The analysis shows that application of the coupled-mode theory in the 
standard form fails to predict the correct scattering properties in case of long crystals. 
In particular, as at short lengths the crystal is reflectionless and invisible when probed 
from one side (according to previous studies [36l Elj), at intermediate lengths the 
crystal remains reflectionless but not invisible. At even longer crystal lengths, both 
unidirectional reflectionless and invisibility properties are broken. 

2. Bragg scattering in T^T-symmetric sinusoidal potentials: general aspects 
and extended coupled-mode theory 

2.1. The model 

Let us consider the stationary Schrodinger equation for a quantum particle in a locally 
periodic and complex potential V(x), which in dimensionless form reads 

H^ = -ft_vix)ij = Ei^ (1) 

where E is the energy of the incident particle and V{x) is the complex scattering 
potential with period A, which is nonvanishing in the interval < x < L. The crystal 
length L is assumed to be an integer multiple of the lattice period A, i.e. L = A^A, 
where N is the number of unit cells in the crystal. As mentioned in the introduction, 
Eq.(l) describes Bragg scattering of matter waves from a complex potential in the non- 
interacting regime, which applies e.g. to a dilute cold atomic beam (see, for instance, 
[29]). In this case, the complex potential arises from the interaction of near resonant 
light with an open two-level system, and it is generally absorptive. In this work we will 
mainly focus our attention to the PT-symmetric sinusoidal potential, assuming 

V{x) = Vo [cos(27rx/A) +iasin(27rx/A)] (2) 

for < X < L, and V{x) = for x < and x > L, where Vq is the lattice amplitude 
and a > measures the strength of the non-Hermitian part of the potential. The 
spectral properties of the T^T-symmetric sinusoidal potential (3) have been investigated 
in Refs.pm [22l l30l [321 [35] . For the infinitely-extended crystal, the energy spectrum 
remains real- valued for a < 1, and breaking of the VT phase is attained at a = = 1 
p!Ol [221 [301 l32l [35] . Here we consider Bragg scattering of incoming waves with momentum 
p close to the Bragg value tt/A, i.e. with energy E = p'^ close to (tt/A)^, and typically 
will assume a modulation Vq of the potential much smaller than the energy E. 
It should be noted that Bragg scattering of optical waves in one-dimensional Bragg 
grating structures, considered in Refs. f36l[37] . is basically analogous to Bragg scattering 
of matter waves in the framework of Eq.(l). In fact, the electric field amplitude S{x) 
of an optical wave at frequency uj that propagates along a dielectric medium with a 
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spatially-dependent relative dielectric constant e{x) = nl[l + Ae{x)]^ where no is the 
refractive index of the lossless medium and Ae{x + A) = Ae(x) accounts for the index 
and gain/loss modulation, satisfies the scalar Helmholtz equation, which can be written 
in the form 

- - EAe{x)£ = E£, (3) 

where we have set E = k'^ and k = noUj/cQ. Note that Eq.(3) formally reduces to the 
Schrodinger equation provided that the following formal substitutions 

S ^ ip , k ^ p (4) 

are made, with a complex scattering potential V{x) related to the modulation of the 
dielectric constant Ae{x) by the simple relation 

V{x) = EAe{x). (5) 

Hence, the only difference between scattering of matter waves in complex optical 
potentials and light waves in complex Bragg gratings is that, in the latter case, the 
complex scattering potential V{x) = EAe{x) in the equivalent Schrodinger equation 
depends on the energy E of the incidence particle. However, for shallow gratings 
Bragg scattering occurs solely for optical fields with frequencies uj very close to the 
Bragg frequency ujb = Co7r/(noA) (see, e.g. [37j), and thus one can safely assume 
V{x) {n / Ay Ae{x) ^ leaving out the dependence of the scattering potential from the 
energy. In the following, we will mainly focus our analysis to the determination of the 
reflection and transmission coeflicients for scattering of matter waves in the framework 
of Eq.(l), however similar results hold mutatis mutandis for reflection and transmission 
of optical waves in Bragg grating structures. 

2.2. Scattering states^ spectrum^ and reflection/transmission coefficients 

Since V{x) = for x < and x > L, the continuous spectrum of the Hamiltonian H is 
the semi-inflnite real axis of energies = > 0, and the corresponding eigenfunctions 
are the scattered states, deflned by the relations 

, . _ J ai ex.p{ipx) + f3i ex.p{—ipx) x < . . 

I Qf2 exp[ip{x — L) + exp[— — L)] x > L 

where p > is the momentum and (o^i, (0^2, /32) are the amplitudes of forward and 
backward propagating waves on the left (x < 0) and on the right {x > L) sides of the 
crystal, respectively. Such amplitudes are related by the algebraic equation (see, for 
instance, [13]) 



where the 2x2 transfer matrix A^(]9) is unimodular, i.e. detA4 = JVi22M.11— = 
1. For a T^T-symmetric potential, the further relation A^22(p) = •^n(]^*) holds. The 
transmission (t) and reflection (r) coeflicients for left (/) and right (r) side incidence are 
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related to the coefficients of the transfer matrix by the usual relations (see, for instance, 

m) 

M22 M22 M22 

Note that the transmission coefficient does not depend on the incidence side like in an 
ordinary crystal, whereas generally one has \r^^^\ ^ i.e. in reffection a complex 

crystal behaves differently for left and right incidence. This is a rather general result of 
wave scattering from a complex potential barrier which was previously discussed e.g. in 
[401 [4T] . If we indicate by Z{p) the fundamental matrix of Eq.(l) from x = to z = L 
which relates the values of and {dip/dx) at the planes x = and x = L, i.e. 

it can be readily shown that the transfer matrix M. can be calculated as 

M{p) = T-\p)Z{p)T{p) (10) 
where we have set 

Up) M . (11) 

y ip —ip J 

From a numerical viewpoint, the fundamental matrix Z{p) can be computed as follows. 
Let us cut the crystal into a sequence of A^o thin slices, of thickness Ax = L/Nq^ and let 
us indicate by Zk{p) the fundamental matrix associated to the propagation at the k-th 
slice {k = 1, 2, A^o), i-e. from Xk = {k — l)Ax to Xk+i = kAx. Then the fundamental 
matrix Z{p) can be calculated as the ordered product 

Z{p) = Zn{p) X Zn-i{p) X ... X Z2{p) X Zi{p) (12) 

If Ax is much smaller than A, the potential V{x) is almost constant in the interval 
(xk^Xk+i)^ and thus Zk{p) can be approximated as 

/ cos(AfcAx) (l/Afc)sin(AfcAx) \ 
^^^^ ~ V -Afcsin(AfcAx) cos(AfcAx) J ^ ^ 

where we have set 



Xk = ^Jp' + Vixk). (14) 

Note that, because of the periodicity of the crystal, one can limit to compute the 
fundamental matrix for the unit cell, Z^^^^^\p)^ i.e. from x = to x = A. The 
fundamental matrix of the crystal is then given by Z{p) = Z^^^^^^^ [p)^ which can 
be computed using the relation [42j Z{p) = Z^^^^^\p)Uni ~ where I is the 

2x2 identity matrix, Un = sin[(A^ + 1)0]/ sin0, and the complex angle 9 is defined by 
cos^ = (l/2)Tr(Z(^^^^)). 

Besides scattering states, the Hamiltonian H can possess bound states belonging to 
the point spectrum, which are determined by the zeros of A422{p) (i-e. the poles of the 
transmission coefficient t) in the half complex plane lm{p) > 0. The zeros of A422{p) 
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in the half complex plane lm{p) < are resonance states. The onset of VT symmetry 
breaking is detected by the appearance of a divergence (for some real value p = po 
of the momentum) in the transmission coefficient as the non-Hermiticity parameter a 
is increased from zero. In fact, for a = the transmission and reflection coeflBcients 
are bounded from above and the zeros of A^22(p) fie in tfie lower fialf complex plane 
lm{p) < 0, i.e. tfiey are resonances. As a in increased, tfie resonances move toward tfie 
real axis lm{p) = 0, until at a critical value a = a resonance crosses tfie real axis, 
say at p = po 7^ (real). Correspondingly, tfie transmission coeflicient t{p) diverges 
dit p = Pq. Just above cTc, a bound state witfi complex energy E = {po + iO^Y tfius 
appears, wfiicfi is tfie signature of VT symmetry breaking. For a sinusoidal crystal of 
finite lengtfi, tfie transmission and reflection coefl&cients at a = 1 are bounded, and 
numerical calculations of tfie transfer matrix Ai (using tfie procedure outlined above) 
indicate tfiat symmetry breaking is attained at a value ac larger thdiU one, witfi ac 1^ 
as L ^ oc. 

2.3. Coupled-mode theory 

For a ratfier general class of PT-symmetric complex potentials V{x) describing Bragg 
scattering in shallow Isittices^ approximate expressions for tfie reflection and transmission 
coeflicients can be derived by an asymptotic analysis of Eq.(l). In tfie optical context, 
sucfi an analysis is generally referred to as tfie coupled-mode tfieory of Bragg scattering, 
wfiicfi is known to provide accurate description of transmission and reflection coeflicients 
for index-modulated sfiallow gratings (see, for instance, [38l |39]). Sucfi an analysis 
applies to Bragg scattering of particles witfi momentum p close to tt/A (flrst-order 
Bragg scattering) provided tfiat tfie particle energy E (tt/A)^ is mucfi larger tfian tfie 
cfiaracteristic modulation deptfi Vq of tfie complex crystal. For tfie sinusoidal crystal 
defined by Eq.(2), tfiis means tfiat tfie parameter a = A^Vq/tt^ sfiould be mucfi smaller 
tfian one. In R.ef . [36l l37] . it was sfiown tfiat application of coupled-mode tfieory to 
tfie sinusoidal PT-symmetric crystal at a = 1 gives tfie following expressions for tfie 
transmission and reflection coeflicients 

t{p) = exp(zj9L), r«(p) = 0, rW(p) = ^^Q^^i^(^-^) exp[i{p+7r/A)L] (15) 

wfiere we fiave set 5 = p — n/A. Sucfi equations clearly indicate tfiat, for left-side 
incidence, tfie crystal appears to be fully invisible, i.e. tfiere are not reflected waves and 
tfie transmitted ones propagate as if tfie crystal were absent [36j . Conversely, for rigfit- 
side incidence as tfie transmitted wave propagates again as if tfie crystal were absent, 
a reflected wave is generated, witfi a reflectance i?^^) = |7-(0|2 ^^j^g^^ grows quadratically 
witfi tfie crystal tfiickness L at Bragg resonance 5 = 0, i.e. for p = tt/A . Sucfi a 
pfiysically relevant befiavior was referred to as unidirectional invisibility in Ref . [36] • As 
it is sfiown in tfie Appendix A and briefly mentioned in Ref. [36] , in tfie framework of tfie 
coupled-mode tfieory unidirectional invisibility is predicted to occur for a ratfier general 
class of complex potentials witfi zero mean, V{x) = X^n/o exp(27rmx/A), provided 
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that the condition $_i = (or similarly $i = 0) is satisfied. Note that the VT- 
symmetric sinusoidal potential (2) at a = 1 belongs to such a general class of complex 
potentials. 

In the next section, we will derive exact expressions for transmission and refiection 
coefficients for the PT-symmetric sinusoidal crystal, and will show that the invisibility 
property of the crystal, as predicted by Eq.(15), occurs solely for short crystals, and 
breaks down for long crystals. Here we discuss the reasons of failure of Eq.(15) 
to predict the correct expressions of transmission and reflection coeflicients in long 
crystals, and propose an extended version of coupled-mode theory to properly describe 
Bragg scattering in complex crystals. The derivation of coupled-mode equations is 
routinely done by means of averaging or multiple-scale asymptotic techniques, which 
is detailed in the Appendix A for the general case of a complex potential V{x) = 
^^^QVnexp{2i7mx/A) with zero mean. Here we give explicit analytical results for the 
sinusoidal T^T-symmetric crystal at a = 1, however as shown in the Appendix A similar 
results are obtained for a more general complex crystal provided that the condition 
$_i = is satisfled. For a small value of a, a solution to Eq.(l) can be searched as a 
power series expansion 

ij{x) = ij^^\x) + a^^^\x) + a^^^\x) + ... (16) 

and multiple spatial scales Xq = x, Xi = ax, ... are introduced to satisfy solvability 
conditions at the various orders in the asymptotic analysis. If the analysis is pushed up 
to the order ^ a, the solution to Eq.(l) inside the crystal for a small value of 5 = p—n/A 
(of order ^ a) can be written as 

ij{x) = (x) + a^^^^ (x) + o{a^) (17) 

where we have set 

ip^^^x) = u{x) exp(i7rx/A) + v{x) exp(— zttx/A) (18) 

a^^^\x) = ^'^'''J'^^ exp(3z7rx/A) (19) 

and where the amplitudes u and v satisfy the coupled-mode equations 

.du VoA . . 

I— = — ou V 20 

rfx 27r ^ ^ 

= Sv (21) 

dx 

From Eqs.(20) and (21), it follows that the amplitudes u and v at the planes x = and 
z = L are related by the relation {u{L)^v{L)Y = /C(]9)(t^(0), 'u(O))^, where the matrix 
lC{p) reads explicitly 

In standard coupled-mode theory [i38l [39], only the leading order term ip^^\x) in the 
expansion (17) is considered for the computation of the transfer matrix and one 
simply has 

M{p) = S{p)lC{p) (23) 



Invisibility in VT -symmetric complex crystals 



8 



where we have set 

^^P>-[ exp(-i7rL/A) ) ' ^^^> 

Using Eq.(23) for the expression of the transfer matrix, one then obtains Eq.(15) for 
the reflection and transmission coefficients. However, for a complex crystal such an 
analysis can fail even if the correction term ailj^'^\ given by Eq.(19), remains smaller 
than To understand why this may happen, let us consider as an example the case 
of left-side incidence. This case requires the absence of incident waves from the right 
side, which formally implies (dijj/dx) = ipip at x = L, where p is the momentum of 
the incident wave from the left side. If the approximation iIj{x) c:^ i/j^^^x) is made, the 
appropriate boundary conditions would be u{0) — 1 and v{L) — 0. The key point is 
that, if the expression (17) of i/j{x) up to the order ^ a is now considered, the boundary 
condition v{L) — does not exactly correspond to the absence of incident waves from the 
right side, just because of the (small) additional contribution to i/j{x) given by aip^^^x) 
[Eq. (19)] . Hence a more accurate procedure should use the boundary conditions u{0) = 1 
and v{L) = e, where e is a small parameter (of order a) to be determined such that 
(diJj/dx) = ipijj X — L. It is obvious that, if the solution to the coupled-mode 
equations (20-21) with the boundary conditions i/(0) = and v{fd) = e (corresponding 
to probing the crystal from the right side with a small amplitude e) would remain small 
uniformly in the interval (0, L), the additional contribution to the solution arising from 
taking e would just introduce a small correction to the solution corresponding to the 
boundary conditions uifi) = 1 and ^'(O) = 0. Hence a small correction to the transmission 
and reflection coeflicients Eq.(15) would be obtained. This case always occurs for short 
enough crystals. However, if the crystal length L is long enough such that the reflection 
^(^) for right-side incidence becomes large (of the order or larger than ^ l/o^), then 
the correction arising from the solution to the coupled-mode equations (20-21) with the 
boundary conditions u{{)) = and ^(0) = e can not be neglected anymore, and thus 
should be accounted for when calculating the reflection and transmission coefficients 
with the appropriate boundary conditions. The crystal length L at which failure of 
Eq.(15) is expected to occur can be estimated by imposing \r^^^ \ '::il/a. Taking for \r^^^\ 
its peak value at 5 = 0, i.e. \r^'^^\ ^ VoAL/(27r), one then expects failure of Eq.(15) for 
L >r^ Lc^ where 

• (25) 



Hence, for L of the order of larger than Lc, Eq.(15) can not be used to calculate the 
transmission and reflection coeflicients of the crystal. A more appropriate procedure 
to compute the transfer matrix, and thus the transmission and reflection coeflicients, 
is to use Eq.(lO), in which the fundamental matrix Z is calculated using for ip{x) the 
expression given by Eqs. (17-19), i.e. including the first-order correction term aijj^^^x). 
Such an extension of the ordinary coupled-mode theory will be referred to as the extended 
coupled-mode theory. 
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3. Bragg scattering in sinusoidal T^T-symmetric crystals: Exact analysis 

Exact expressions for the reflection and transmission coeflicients can be derived for the 
sinusoidal potential (2) at a = 1 in terms of modifled Bessel functions of Arts kind. 
In fact, after the change of variable y = (Aa/T^/tt) exp(i7rx/A), Eq.(l) reduces to the 
Bessel equation [35] 

+,)^ = (26) 

where we have set 

, = '-^^ (27) 

TT 

Note that Bragg scattering for particles with momentum p close to tt/A corresponds to 
q ^ 1. For g 7^ 1, two linearly independent solutions to Eq.(26) are Iq{y) and I-q{y)^ 
where Iq{y) is the modifled Bessel function / of flrst kind [43]. Hence in the interval 
< X < L two linearly independent solutions to Eq.(l) are given by 

$i(x) = /^(Aexp(i7rx/A)) , $2(^) = /-g(Aexp(i7rx/A)) (28) 

where we have set 



a. (29) 

TT 

Using Eq.(28), one can construct the fundamental matrix Z[x) of Eq.(l) from x = to 
X ^ L as 

^^'^ - [ $;(L) $;(L) ) [ $;(o) m ) ^''^ 

where the apex denotes the derivative with respect to x. The transfer matrix A4{p) 
is flnally obtained after substitution of Eq.(30) into Eq.(lO). After some lengthy 
calculations, which are briefly detailed in the Appendix B, the following expressions 
for the transfer matrix coeflicients are obtained: 

. Asin(]9L) 

2p sin(7rg) 
Asm{pL) 



Mn{p) = cos(pL) + 1-^^^^ {p'QiQ2 - VoD.D^) (31) 



Mi2{p) = -i 



M2iip) = i 



2psin(7rg) 
A sm{pL) 
2psm{7rq) 



V0D1D2 + p'QiQ2 + pJVo {D1Q2 + D2Q1) 



V0D1D2 + p'QiQ2 - P^Vo {D1Q2 + D2Q1) 



(32) 
(33) 



M22{p) = cos(pL) - I . [ (p'QiQ2 - VoD,D2) (34) 

2psm(7rg) ^ ^ 

where we have set 

= /,(A) , Q2 = /-,(A) , D, = /;(A) , D2 = / ,(A) (35) 

and where q and A are deflned by Eqs.(27) and (29), respectively. The reflection and 
transmission coeflicients r^^'^\p) and t{p) are then obtained after substitution of Eqs.(31- 
34) into Eq.(8). In particular, for the transmission coeflicient t{p) one obtains explicitly 

tip) — —T — — — — - (36) 

co8[pL) — iF{p) 8m{pL) 
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where we have set 

FiP) = ^ , (p'QiQ2 - VoD.D^) . (37) 

Note that Eq.(36) would reduce to the first of the Eq.(15), corresponding to 
unidirectional crystal invisibility, if the function F{p) were replaced by 1. 
The previous equations (31-37) have been derived for Bragg scattering of matter waves 
in the framework of the Schrodinger equation (1), however similar relations hold for 
Bragg scattering of light waves in a complex Bragg grating structure, governed by the 
similar equation (3). Specifically, in view of Eqs.(4) and (5), for a grating structure 
with a sinusoidal PT-symmetric modulation of the complex relative dielectric constant 
Ae{x) = $ exp(2i7rx/A) and for incident light waves with frequency cj, the analytical 
expressions given above still hold, provided that the following formal substitutions 
p (jono/co and Vq (nocj/co)^^ — (7r/A)^$ are made, where no is the refractive 
index of the lossless dielectric medium and Cq the speed of light in vacuum. 



4. Unidirectional crystal invisibility 

Let us know discuss the unidirectional invisibility of the sinusoidal PT-symmetric 
crystal on the basis of the exact scattering results presented in the previous section. 
To study the exact behavior of t{p) as given by Eq.(36) and breakdown of crystal 
transparency as the number of cells is increased, let us measure the length x 
in units of A/tt, i.e. let us set without loss of generality A = tt. With such a 
scaling, one has a = Vq^ q = p^ L = Nn and A = y^. Using the identity 
/^(A) = Vi(A) - {p/A)Ip{A) = Vi(A) + {p/A)Ip{A) for the derivative of modified 
Bessel functions j43], one can write 

^^^^ " cos(A^7rj9) - iF{p) 8m{N7rp) ^^^^ 



with 

TT 



Fip) 



(39) 



2p sin(7r]9) 

and where the Bessel functions are calculated at ^/V^ . Note that t{p) depends on 
two parameters solely: the number of crystal cells and the potential amplitude 
Vq. Note also that Eq.(38) is valid regardless of the smallness of Vq and far from the 
Bragg resonance condition p = 1 as well. Similar expressions can be derived for the 
refiection coefficients r^^'^\p) in terms of modified Bessel functions. In the computation 
of the reffection and transmission coefficients, the modified Bessel functions have been 
calculated using a fast and highly accurate routine, discussed in Ref. [44j. The accuracy 
of our procedure has been tested by checking the agreement of the results obtained 
from the exact analytical prediction [Eqs. (38) and (39)] and from the full numerical 
procedure outlined in Sec. 2. 2 [Eqs. (10-14)]. 

Figures 1 and 2 show the behaviors of spectral transmittance T{p) = |t(p)P and 
refiectances R^^'^\p) = \r^^'^\p)\'^ for left and right side incidence, as calculated by the 
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Figure 1. Behavior of (a) transmittance, (b) reflectance for left-side incidence, (c) 
normalized phase time (in transmission), and (d) reflectance for right-side incidence 
versus momentum p of the incident particle in a sinusoidal PT-symmetric crystal for 
parameter values A = tt, = 0.02 and for a number of cells = 50. The curves are 
obtained by using the exact expressions Eqs. (31-34) for the transfer matrix coefficients. 



exact analysis (solid curves), for increasing values of the number of cells and for 
Vo = 0.02. In the figures, the behavior of the normalized phase time of transmitted 
waves, defined by rt{p) = {1/ L){d(j)t/ dp)^ is also depicted, where (j)t{p) is the phase 
of t{p). Physically, Tt{p) represents the traversal time of a narrow wave packet, with 
central momentum across the crystal, normalized to the transit time in vacuum (i.e. 
in the absence of the crystal). For a relatively small numbers of cells, as in Fig.l, 
unidirectional invisibility is observed, and reflection for left-side crystal incidence is 
extremely small, according to the coupled-mode theory of Sec. 2. 3 and the results of 
Refs.[36i l37] . However, as the number of cells is increased to become comparable or 
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Figure 2. Same as Fig.l, but for N = 2000. In the figures, solid lines refer to the 
exact scattering results obtained from Eqs. (31-34), whereas the dashed curves are the 
predictions based on the extended coupled-mode theory. 



larger than Nc = Lc/A^ given by [see Eq.(25)] 

iVe ^ ^ = (40) 

the invisibility regime breaks down near p = 1^ with the appearance of oscillations 
of the transmittance and phase time, as one can clearly see in Figs. 2(a) and (c). In 
the figures, the predictions of spectral transmittance and reflectance computed by the 
extended coupled-mode theory, discussed at the end of Sec. 2. 3, are also depicted by the 
dotted curves. Note that, as the crystal is not anymore transparent, the reflectance 
for left-side incidence remains extremely small [see Fig. 2(b)]. In such a regime the 
crystal is not invisible, however it is still unidirectional reflectionless. As the number 
of cells is further increased, the transmittance grows in a narrow interval near the 
Bragg condition p 1^ slS well the the reflectance for right-side incidence. As the 
transmittance becomes large enough, such that A^22 becomes smaller and of the same 
order of magnitude than A^2i, a very narrow resonance peak appears in the reflectance 
spectrum for left-side incidence. This is shown in Fig. 3, in which solid and dotted 
curves refer to the exact results and to the approximate ones based on the extended 
coupled-mode theory, respectively. In such a regime, both unidirectional invisibility and 
reflectionless properties of the complex crystal are thus broken. To estimate the number 
of cells N'^ > Nc at which such a second transition occurs, we can apply the extended 
coupled-mode theory, discussed in Sec. 2. 3, to calculate an approximate expression of the 
reflection coefficient r^^^ for left-side incidence. Using Eqs. (17-21) for an approximate 
expression of '0(x) and applying the appropriate boundary conditions, corresponding to 
left-side incidence, at exact Bragg resonance (p = tt/A) an approximate expression for 
|r(^)| can be derived, which reads explicitly r^^^ ^ (7r/64)a^(L/A). The critical number 
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of cells A^^ at which the crystal is not anymore reflectionless for left-side incidence can 
be estimated by letting |r(^^| ^ 1, which yields 
64 

K - —r (41) 




Figure 3. Transmittance (left panels) and reflectance for left-side incidence (right 
panels) in a sinusoidal "PT- symmetric crystal for Vq = 0.02, A = tt and for increasing 
number of crystal cells: (a) N = 10000, (b), 50000, and (c) TV = 1600000. 



5. Conclusion 



In this work Bragg scattering in sinusoidal PT-symmetric complex crystals of finite 
thickness has been theoretically investigated, and exact analytical expressions for 
reflection and transmission coeflicients have been derived in terms of modifled Bessel 
functions. The analytical results indicate that unidirectional invisibility, recently 
predicted for such crystals by coupled-mode theory |[36l [37], breaks down for crystals 
containing a large number of unit cells. In particular, for a given modulation depth 
in a shallow sinusoidal potential, three regimes have been found as the crystal length 
is increased. At short lengths the crystal is reflectionless and invisible when probed 
from one side (unidirectional invisibility), according to standard coupled-mode theory. 
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As the numbers of cells is increased, absence of reflection for one side incidence is 
still observed, however the crystal is no more invisible because large oscillations in the 
transmittance and in the transmission phase time appear near the Bragg resonance. 
For still thicker crystals, both unidirectional reflect ionless and invisibility properties are 
broken. An extension of coupled-mode theory has been proposed to properly modelling 
the scattering properties of complex crystals. 



Appendix A. Derivation of coupled-mode equations 

In this Appendix we briefly derive coupled-mode equations describing flrst-order Bragg 
scattering in a complex crystal with a shallow lattice in the framework of the Schrodinger 
equation (1). Let us consider Bragg scattering of a particle with momentum p close to 
the Bragg value tt/A, i.e. with energy E = p^ close to (tt/A)^, and let us assume for the 
complex scattering potential V{x) a rather general profile with zero mean, given by the 
Fourier expansion 

V{x) = ^ $n exp(2i7rnx/A) (A.l) 

for < X < L. The shallow lattice approximation implies that the Fourier amplitudes 
$^ are much smaller than E. Note that the PT-symmetric sinusoidal potential (2) 
is obtained as a special case of Eq.(A.l) after setting $i = (Vo/2)(l + a), $_i = 
(Vo/2)(l — cr) and = for n 7^ ±1. To develop a perturbative analysis of Eq.(l) in 
the shallow grating approximation V{x) and for p ^ tt/A, it is worth introducing a 
parameter a that measures the smallness of the various terms entering in the equations 
and rewriting Eq.(l) in the following form, suited for an asymptotic analysis 

where we have set 



+ (x) V' = -«[n^)V' + W^V'], (A.2) 



The problem is to construct an asymptotic approximation of the perturbed solution 
= iIj{x] a) to Eq.(A.2) as a ^ 0. Therefore, we seek a perturbation expansion of '0 in 
the form of a power series in a 

a) = i;^^\x) + ^^^^(x) + a^i;^'^\x) + ... (A.4) 

and introduce multiple scales for space x 

Xq = X , Xi = ax , 5X2 = a^x , .... (A. 5) 

which are needed to satisfy the solvability conditions in the asymptotic expansion at 
various orders. Substitution of the Ansatz (A.4) into Eq.(A.2) and using the derivative 
rule 

(f ^ d d of 9^ ^ d d \ , , 

^ = axf + '"5XoaX^ + " [m^^dXodxJ^''' ^^-'^ 
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yields a hierarchy of equations for successive corrections to '0, which are obtained after 
collecting the terms of the same order in a in the equation so obtained. At leading order 
^ one has 

whose general solution is given by 

(^0,^1, ^2,...) = u{X^,X2, ...)exp(Mo/A) + 

+ v{Xi, X2, ...) exp(-i7rXo/A) (A.8) 

where the amphtudes u and v may vary over the slow spatial scales Xi, X2, ... At order 
~ a, one obtains 



dX^ \AJ 
where the forcing term g^^"^ is given by 

/92 / (O) 

5W(Xo) = - [V{Xo) + W] - 2^^^. (A.IO) 

The solvability condition for Eq.(A9) requires that the forcing term g^^\Xo) does not 
contain terms oscillating like ^ exp(±i7rXo/A). After substitution of Eqs.(A.l) and 
(A.8) into Eq.(A.lO) and letting equal to zero the coefficients of the terms oscillating 
like ^ exp(±i7rXo/A) in the expression so obtained, the following equations for the 
evolution of the amplitudes u and v on the slow spatial scale Xi are then obtained 

2%-—- ^ -Wu- ^iv (A.ll) 
koXi 

TT dv 

2i--—^Wv + ^.iu (A.12) 
A 0X1 

and the solution at order ~ o; is given by 

$nexp[^7^(2n + l)Xo/A] 

„ ■i»„expM2n-l)Xo/A| 
+ ^"„|5, (2n-l)'-l 

The evolution equations of the envelopes u and v at longer spatial scales X2, X3, ... 
are obtained similarly as solvability conditions at orders a^, a^, .... in the asymptotic 
expansion. The evolution of the amplitudes u and v in the physical spatial variable x 

are then given by du/dx = adx^u + a^d\^u + .... and dv/dx = adx^v + a^d\^v + If 

we limit our analysis to the order ^ a, after setting a = 1 from Eqs.(A.ll) and (A.12) 
and using Eq.(A.3) one finally obtains the following coupled-mode equations for the 
envelopes u and v 

.du^ -^iz-pi?; (A.14) 
dx 

i-— — 6u + P2U (A. 15) 

dx 
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where we have set 
6 = p 



Pi 



A$i 



P2 



(A.16) 



In the framework of the standard coupled-mode theory [36l [37], unidirectional crystal 
invisibility is attained whenever the evolution of either one of the two amplitudes u oi v 
is decoupled from the other one, i.e. for either $_i = or $i = 0. In particular, for the 
PT-symmetric crystal [Eq.(2)] at a = 1, one has $i = Vq and $_i = 0, which ensures 
unidirectional invisibility. In this case, the coupled-mode equations (A. 14) and (A. 15), 
as well as the first-order correction t/;^-^^ as given by Eq.(A.13), reduce to Eqs.(19), (20) 
and (21) given in the text. 



Appendix B. Derivation of the transfer matrix 

The exact expression of the transfer matrix A4{p) is obtained using Eq.(lO), where the 
fundamental matrix Z{p) is calculated in terms of modified Bessel functions according 
to Eq.(30). A simplified form of Z{p) can be obtained after observing that, owing to 
the analytic continuation of the Bessel / function in the complex plane [43j , one has 

$i(x+L) = $i(x) exip{i7rqL/A), ^2{x+L) = $2(^) exip{-i7rqL/A) (B.l) 

and thus $i,2(^) = $1,2(0) exp(±i7rgL/A) and $i^2(^) = ^[^2{^) ew{^i^QL/A). Hence 
one can write 

/ $2(L) \^( $1(0) $2(0) A / exp(i7rgL/A) 

V $;(L) ) [ $;(0) $;(0) ) [ O exp(-.7rgL/A) 

Moreover, taking into account the property of the Wronskian of modified Bessel 
functions [43j, one has 

$i(0) $2(0) ^ ,7rA /,(A) /_,(A) 

$;(o) $;(o) 'a /;(a) 



-i- sin(g7r) 



(B.3) 



and hence 



$1(0) $2(0) 
$;(o) $;(o) 



A 



= I- 



$;(o) 



-$2(0) 
$1(0) 



(B.4) 



'2sin(g7r) V -$i(0) 

Substitution of Eqs.(B.2) and (B.4) into Eq.(30) yields for the fundamental matrix Z{p) 
the following simplified expression 

= ^ ( *2(0) \ ( exp{i7rqL/A) \ 

2sin(g7r) 1^ ^[{0) $^(0) / V exp{-iTrqL/A) J 

$'2(0) -$2(0) \ 
-$;(0) $i(0) ) 

where $i(0) = Qi, $2(0) = Q2, $i(0) = iVVoD^, $2(0) = iV^)D2 and Qi,2, A,2 are 
defined by Eq.(35) given in the text. Substitution of Eq.(B.5) into Eq.(lO) and using the 
expression of the matrix T{p) given by Eq.(ll), after some lengthy but straightforward 
calculations one finally obtains for the coefficients of the transfer matrix A4{p) the 
expressions given by Eqs. (31-34). 



X 



(B.5) 
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